Structural compliance, misfit strain and stripe nanostructures in cuprate 

superconductors 
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Structural compliance is the ability of a crystal structure to accommodate variations in local atomic 
bond-lengths without incurring large strain energies. We show that the structural compliance of 
cuprates is relatively small, so that short, highly doped, Cu-O-Cu bonds in stripes are subject 
to a tensile misfit strain. We develop a model to describe the effect of misfit strain on charge 
ordering in the copper oxygen planes of oxide materials and illustrate some of the low energy stripe 
nanostructures that can result. 
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The existence of charge stripes in high T c superconduc- 
tors [p] and their relevance to the high-T c property [^]]|] 
are hotly debated subjects at present. The tendency to- 
wards charge phase separation and stripe formation has 
been seen in models of strongly correlated electron sys- 
tems |^|-|| . Charge stripes have also been seen directly in 
strongly correlated systems such as manganites |tJ and 
nickelates ||. In the cuprates, convincing direct evidence 
for stripes comes with the observation of charge-order 
superlattice peaks in Lai.g_ 2; Ndo.4Sr a ;Cu04, though the 
observation of static long-range ordered stripes correlates 
with poor or nonexistent superconductivity |l]|| . The in- 
terest has shifted to dynamic fluctuating short-range or- 
dered stripes. Indirect experimental evidence exists from 
inelastic neutron scattering supporting their presence in 
La 2 _ 3; (Sr,Ba) 2; Cu04 HI and YBa 2 Cu 3 6+( 5 012). The 
observation of a pseudo gap in the electronic density of 
states above T c is even more widespread fi~3|| and can 
be interpreted as originating from dynamic stripes. Of 
critical importance is an understanding of what prevents 
stripes ordering over long-range and becoming static and 
insulating. Here we stress the importance of lattice strain 
to this phenomenon. 

Strong coupling of doped charges to the lattice is 
clearly evident in the cuprates [|l ". 
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l^Jl5[| . A direct way 
in which the charge couples to the lattice is through the 
Cu-0 bond length which depends on charge filling of the 
Cu-0 er* covalent bond |l(|. Doping holes into the Cu02 
plane reduces charge density in this bond which stabi- 
lizes and shortens the bond. This is seen experimen- 
tally as a shortening of the Cu-0 bond with increasing 
doping [ pT[ and has profound implications when charge- 
stripes are present. The presence of charge stripes im- 
plies that regions of the Cu0 2 plane that are heavily 
doped, and therefore have short bonds, coexist with un- 
doped regions with long bonds. Experimental evidence 
for this coexistence has been presented in superconduct- 
ing La2- K (Sr,Ba)a;Cu04 samples Q. We show that this 
results in a misfit strain that breaks up the stripes into 
a domain microstructure. The characteristic length-scale 
of the domains depends on the misfit strain and can range 



from nanometer to long range. The importance of elas- 
tic compliance and interface strain to texture formation 
has been discussed in related ferroelastic materials |nj 
but this is the first time it is discussed in the context of 
breaking up stripes and forming stripe nanostructures. 

A key property of perovskite structures is their capac- 
ity to accommodate atoms of various sizes. They do this 
via a variety of buckling and tilt modes. These modes 
have relatively small strain energy costs, and are fre- 
quently called soft or floppy modes [ ^6|j20|pl[ | . In man- 
ganites, nickelates and cuprates, the key soft mode is a 
buckle in the Cu-O-Cu bond and doping typically relaxes 
this buckle in these materials. This low energy accom- 
modation affects the stability of the structure against 
polaron formation |22|]. We argue that in cuprates, the 
short Cu-O-Cu bond is not fully accommodated by relax- 
ation of the buckle present in the undoped copper-oxygen 
bond. This leads to tensile misfit strains in the copper 
oxygen planes of cuprates. In the presence of stripes the 
tensile misfit strains can lead to a variety of domain mi- 
crostructures. We develop a model to describe these mi- 
crostructures and present some illustrative results using 
the model. 

The ability of the structure to accommodate local bond 
shortening we call "structural compliance" 5 S . Where 
this accommodation comes from octahedral tilting it is 
defined as 5 S — — (\ _ cos a) where Tb and Tf 

are the lengths of the Cu-0 bond in the "buckled" and 
"flat" configurations respectively and a is the angle the 
Cu-0 bond makes with the plane (Fig. [l]). In the reduced 
I4/mmm unit cell of the La2- £C (Sr,Ba) 2; Cu04 system r/ is 
simply a/2 where a is the lattice parameter. In general 
a is small and S s w a 2 /2. The structural compliance 
therefore depends very sensitively on the buckling angle. 
In the case of low-temperature orthorhombic (LTO) p7[ , 
and tetragonal (LTT) 23 1, octahedral tilts observed in 
La 2 _ a; (Sr,Ba) a; Cu04, a = 9 LTT = 9 LTO /V^ where 6 LT t 
and 6 lto are the octahedral tilt angles in the LTT and 
LTO phases respectively. We note that for the same tilt 
angle the LTT symmetry tilts have twice the structural 
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compliance. 

However the La2- s (Sr,Ba)a;Cu04 system has a rela- 
tively small structural compliance, as do the other high 
T c cuprates. On going from undoped to heavily doped 
(x = 0.20) in the La2- 2: (Sr,Ba) 2 ;Cu04 system at 10K, the 
Cu-0 bond shortens from r b = 1.904 A to r b = 1.882 A, 
i.e., by ~ 0.022 A frL7| . For illustration purposes, we 
take this bond length difference to be an estimate of the 
difference in length of the short and long Cu-0 bonds 
in the Cu-0 planes. In order for structural compliance 
to accommodate this bond length mismatch we need 
5 S = 0.022/1.89 = 0.012, which requires a tilt angle of 
8.9° in the LTT symmetry and of 12.6° in the LTO sys- 
tem. However in the La2- a; (Sr,Ba) 2 ;Cu04 system tilts 
have LTO symmetry with a tilt angle of ~ 3.5° resulting 
in a structural compliance of only S s ~ 0.001. There is 
therefore insufficient structural compliance to accommo- 
date the short Cu-0 bonds in these materials, resulting in 
a misfit strain if charges localize, for example, in stripes. 
Note that these numbers represent rough estimates and 
the precise value for the crossover will depend on the 
charge density and degree of derealization of the charge 
in the stripe as well as details of whether the stripes are 
site or bond centered. Nevertheless, it is evident that in 
cuprates, stripe formation will lead to misfit strain, but 
that materials with moderate tilts could accommodate 
short bonds with little or no misfit strain. 

We now consider what happens to the microstructure 
of the charge stripes in the case where a misfit strain 
is present. We expect/show that the presence of strain 
results in a break-up of the stripes into short segments. 
Lattice strain therefore prevents the stripes from order- 
ing over long range; a situation that is a prerequisite for 
most stripe theories of high temperature superconduc- 
tivity P,p ,p4 - 28 1 . This has implications for supercon- 
ductivity by producing an electronic microstructure ]29| ] 
of domains of broken stripes that results in a high den- 
sity of topological stripe defects p5pc| , |3l| ] . The ends of 
the stripes are in a state of high stress making them es- 
pecially susceptible to charge fluctuations. 

We assume that neighboring doped Cu sites in a stripe 
have a separation, l that is shorter than, but approxi- 
mately equal to, the average separation, a, dictated by 
the crystal structure and the average doping level. In this 
case, as the stripes form and increase in length strain en- 
ergy will build up. The basic competition is between the 
length scale of the periodic potential a due to the lat- 
tice structure and the natural length of the short bonds 
due to near neighbor doped sites, as illustrated in Fig. [|. 
The breakup length scale of the stripes is approximated 
by comparing the strain energy gain (~ k nn N(a — lo) 2 /2) 
which occurs when a strained stripe of length L = Na 
breaks up, to the formation energy per bond of an un- 
strained stripe (— J nn ). This yields (dropping constant 
prefactors) 



The key parameters setting the scale of the energies in 
Eq. (1) are k nn , which is proportional to the stiffness of 
the Cu-0 bonds and J nn the energy gain (per bond) in 
forming stripes. The spring constant k nn ~ Vq/o 2 , where 
Vb is of the order of an ionic bond energy (~ 5 eV). 
J nn is of the order of the stripe formation tempera- 
ture, e.g., the temperature where the pseudogap appears 
~ 400K ~ AOmeV. As discussed above, the misfit 
strain for cuprates is of order (a — lo)/a ~ 0.01. Us- 
ing these numbers in Eq. (1) indicates that for cuprates, 
the characteristic length of stripe domains, N c , should 
be of order 100 lattice spacings. The stripe length di- 
verges quadratically (see Eq. (1)) as the misfit strain 
approaches zero leading to long range stripes in oxides 
with no misfit strain, as appears to be the case, for ex- 
ample, in La2-ir-yNdj,Sr :E Cu04. 

We now develop a model to test whether the simple 
prediction (1) holds when there are many stripes and 
to study the stripe nanostructures produced once strain 
destabilizes long-range stripes. This model is as simple 
as possible while including the key aspects of doping, 
stripes and misfit strain. Note that this model consid- 
ers only stripe breakup and not stripe formation. We 
assume magnetic effects and long-range Coulomb inter- 
actions play a key role in stripe formation [||j3^j3^], but 
that misfit strain leads to their breakup. Quantum fluc- 
tuations may also contribute to stripe breakup, however 
we consider here only misfit strain which we assume is 
the dominant factor. The model is defined on a square 
lattice where each site is assigned an occupancy number 
Hi where rij = 1 indicates that the site is doped with a 
hole and rij = indicates that it is undoped. 



H= - 2^ij JnnniUj +2^ij JnnnniUj + 

+ Knikj ~ lo) 2 nin 3 + J^i riiV{fi). 



(2) 



(1) 



In the spirit of the stripes seen in Lai.6-xNdo.4Sr :E Cu04 
where a hole resides on every second site along the stripe, 
we can simplify the computation by considering a sub- 
lattice of every second Cu site and include an attractive 
nearest neighbor and repulsive next nearest neighbor in- 
teraction, J nn and Jnnn respectively. These will result in 
stripe formation provided J nn > 0, J nnn > 0. The third 
term in the Hamiltonian has the chemical potential, /x, 
and allows us to vary the doping in the model to some 
occupation density, p. Note that our sublattice model 
only allows doping at every second copper site along the 
x and y directions so that p is related to the doping frac- 
tion p in the cuprates via, (p) « 4p. The last two terms 
in Eq. ^ are the strain terms. The first term contains 
the misfit strain which occurs when two nearest neighbor 
sites on the sublattice are occupied (this corresponds to 
next-nearest neighbor sites in the cuprates). The param- 
eter Iq is the natural length of this bond, i.e., the length 
that bond would take in the absence of stress. The sec- 
ond term V(fi) is the periodic potential imposed on the 
copper oxide layers by the crystal structure. In our model 
we take this potential to be modulated along both of the 
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in-plane axes of the CuC>2 planes with a wavelength a, 
and to have a amplitude Vq. These terms in the Hamil- 
tonian result in a finite stress energy unless Iq = a as can 
be seen in the schematic in Fig. |2[ The strain energy is 
long range and has to be relaxed self consistently using 
a gradient or Newton's method. 

It is difficult to find the ground state of the strained 
lattice gas (2), however it is easy to demonstrate that 
long range stripes are unstable in the presence of mis- 
fit strain and to test the general form (1). To illustrate 
this, consider the three stripe nanostructures of Fig. ||. 
From the lattice gas model (2), it is evident that the long 
range stripe has the lowest energy, provided there is no 
misfit strain. However in the presence of misfit strain, 
the lower two nanostructures of Fig. |^ have lower strain 
energy, though they do incur an energy cost proportional 
to J nn due to breaking up the long range stripes. The 
strain energy in the long stripes is simply k n n(la — c) 2 /2 
per bond. For a strain of 0.01 and Vq = 1 = k nn (in lat- 
tice units), this leads to a strain energy of 0.0152 for the 
structure of the top panel in Fig. ||. In contrast the inter- 
leaved stripe and weave structures (bottom two panels in 
Fig. |J) have strain energies of 0.0100 and 0.00765 respec- 
tively. These calculations were carried out using simple 
gradient descent on the unit cells of Fig. ||, with the peri- 
odic cell shape fixed, and the boundary atoms fixed and 
excluded from the total energy calculation. Tests using 
free and periodic boundaries and other nanostructures 
indicate that the state of lowest energy does depend on a 
variety of factors, including doping, boundary conditions 
and cell size. However in all cases, the strain energy of 
the long range stripes is of order twice that of the stripe 
domain states and this is not sensitive to the modeling 
details. It is thus reasonable to find the breakup length 
scale of the long range stripe array (top panel of Fig. ||) 
by comparing the strain energy of long-range stripes with 
the nearest neighbor energy cost incurred in breaking the 
stripes. This leads to N c k nn (lo — a) 2 m J nn , which is 
consistent with the simple estimate (1). In general, the 
ground state structures of discrete strained systems are 
extremely rich [|34j-|38| and we expect a similar richness 
in the ground structure of the Hamiltonian (2) as a func- 
tion of doping and temperature, as will be described else- 
where. However the main point emphasized here is that 
misfit strain destabilizes long range stripes in cuprates at 
a length scale well approximated by (1). 

Now we briefly comment on a number of experimen- 
tal observations that can be explained by strain induced 
stripe nanostructures. The most important point is that 
the presence of strain will prevent the stripes from or- 
dering over long range resulting in broken stripes even 
at moderately high dopings. Broken stripes will have 
significant quantum fluctuations. More structural com- 
pliance in the form of enlarged octahedral tilt distortions 
will reduce the misfit strain and cause the stripes to grow 
in length, presumably slowing their dynamics. Systems 
where the Cu02 plane is in tension, or nearly in tension, 
such as the Hg and Tl systems will have faster stripe dy- 



namics and therefore better superconducting properties 
and higher T c 's as observed. Similarly, the increase in T c 
with decreasing tilt distortions (3^^T|| have an explana- 
tion here. The strained stripes will also tend to stabilize 
the LTT phase over the LTO phase because this results in 
greater structural compliance, as we discussed, and there- 
fore a decrease in the misfit strain. Stripe microstruc- 
tures such as those shown in the third panel of Fig. || 
will result in local LTT symmetry tilts but averaging 
over the two 90° rotated LTT variants will lead to effec- 
tive LTO behavior, as observed in La2_ K Ba K Cu04 pl| . 
Also, severely misfitting ions will give rise to large ampli- 
tude tilts locally, resulting in locally varying structural 
compliance, which will tend to pin the stripes. This ex- 
plains the observation that superconductivity can be sup- 
pressed by increasing the mean-square dopant-ion size at 
fixed average dopant-ion size [0. Finally, we note that 
ferroelastic effect s Ba| and the observation of tweed mi- 
crostructures [ |19[|35||1 4 1 , as well as the sensitivity of T c 
on dopant ordering in Yf^CuaOe+a and La2-a;Cu04+5, 
might have an explanation in strained stripe induced mi- 
crostructures. 

In summary, we point out that doping shortens the 
Cu-0 bond length and that inhomogeneous doping in 
the form of stripes results in a tensile misfit strain. This 
breaks up the stripes and results in, presumably fluctu- 
ating, domains of broken stripes where the domain size 
depends on the misfit strain. Structural compliance in 
the form of buckling modes can reduce the misfit strain 
and result in long range ordered and static or quasi-static 
stripes reducing or destroying superconductivity. 
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FIG. 1. The angle a which describes the buckle in the cop- 
per (solid circle) - oxygen (open circle) - copper bonding. 
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FIG. 2. Illustration of the effect of tensile misfit strain on 
structure. The periodic potential has wavelength a, the lattice 
spacing set by the crystal, while the natural length of the 
shorter doped Cu-Cu bonds is lo < a. Due to the tensile 
misfit, the Cu atoms are unable to remain at the minima of 
the periodic potential and, at the same time, to conform to 
the natural length Iq. 
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FIG. 3. Stripe nanostructures in cuprates. In systems with 
tensile misfit, the static stripe (top panel) has a high strain 
energy. Interleaved stripes (middle panel) and the weave mi- 
crostructure (bottom panel) reduce this strain energy as they 
can accommodate strain relaxation. One unit cell of each mi- 
crostructure is shown. Only copper atoms are shown. Con- 
centric circles indicate doped sites, while sites indicated by 
crosses are excluded from the strain relaxation. For illustra- 
tion purposes the natural length of bonds between doped sites 
are 10% shorter than bonds between undoped sites. 
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